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Abstract 

Resistive steady states in toroidal magnetohydrodynamics (MHD), where 
Ohm's law must be taken into account, differ considerably from ideal ones. 
Only for special (and probably unphysical) resistivity profiles can the Lorentz 
force, in the static force-balance equation, be expressed as the gradient of a 
scalar and thus cancel the gradient of a scalar pressure. In general, the Lorentz 
force has a curl directed so as to generate toroidal vorticity. Here, we calculate, 
for a collisional, highly viscous magnetofluid, the flows that are required for an 
axisymmetric toroidal steady state, assuming uniform scalar resistivity and 
viscosity. The flows originate from paired toroidal vortices (in what might 
be called a "double smoke ring" configuration), and are thought likely to be 
ubiquitous in the interior of toroidally driven magnetofluids of this type. The 
existence of such vortices is conjectured to characterize magnetofluids beyond 
the high-viscosity limit in which they are readily calculable. 
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I. INTRODUCTION 



If static solutions of the MHD equations in toroidal geometry are forced to obey Ohm's 
law as well as force balance, most of them disappear. The reason is that the current den- 
sities resulting from inductively generated, steady electric fields, give the term j x B in the 
equation of motion (where j is the electric current density and B is the magnetic field) a 
finite curl, and thus j X B cannot be balanced by the gradient of a scalar pressure. Q The only 
exceptions, for finite, uniform transport coefficients and incompressible MHD, involve resis- 
tivity profiles which are not uniform on a magnetic flux surface, and thus may be thought 
unrealizeable.i Earlier, it was speculatedB that if toroidal MHD steady states exist for resis- 
tive magnetofluids, they will probably involve flows (vector velocity fields). Here, we wish to 
calculate and describe such velocity fields for the case of large plasma viscosity (low viscous 
Lundquist number). We stress that the effect we are describing is a consequence of toroidal 
geometry, and is not an issue in the "straight cylinder" approximation. Vortices similar 
to the ones found here can appear in straight-cylinder MHD computations but only above 
instability thresholds.^ In the toroidal case, the vortices are an integral part of the force 
balance in the steady state all the time, and are not uniquely connected with instabilities. 

In Sec. II, the governing MHD equations are written out, and an approximate method 
of solution valid in the limit of large viscosity (low viscous Lundquist number) is discussed. 
With our approach, the velocity fields necessary to maintain axisymmetric, toroidal, resis- 
tive steady states are determined. In Sec. Ill, numerical implementation of the procedure 
is described, and vorticity contours and streamlines resulting from the calculations are pre- 
sented. We believe the results may suggest a simultaneous occurrence of current circuits 
and vortex rings in a finitely electrically conducting fluid that goes well beyond this single 
tractable example. Sec. IV is a summary and suggestion for further investigations. The Ap- 
pendix describes, in terms of inequalities among dimensionless Reynolds-like numbers, the 
approximations made and also specifies the geometry in detail. Purely for computational 
convenience, we specialize to a toroid with a rectangular cross section (see the Appendix 
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and Fig. 1). However, we believe our results to apply to toroids with more general boundary 
shapes. 

II. A HIGH- VISCOSITY CALCULATION 

The MHD equation of motion for a uniform-density, incompressible magnetofluid is@ 

— + v • Vv = jxB - Vp + zA7 2 v, (1) 
at 

where v is the fluid velocity, p is the pressure, and v is the kinematic viscosity. We work in 
standard "Alfvenic" dimensionless units, so that in fact v is the reciprocal of a Reynolds- 
like number; specifically, v~ x is the viscous Lundquist number, given in terms of quantities 
expressed in cgs units by 

v = C a L/D, (2) 

where C a is an Alfven speed based on the mean poloidal magnetic field, L is a characteristic 
length scale which can be taken to be a toroidal minor radius, and v is the laboratory 
kinematic viscosity, expressed in cm 2 / s. We assume that the magnetofluid is incompressible, 
and that the viscosity and electrical conductivity are spatially uniform scalars. These are 
all significantly restrictive assumptions, but they will be seen to lead to a tractable problem 
in otherwise uncharted territory. 

The most severe assumption we will make is that of small viscous Lundquist number. 
For a collision-dominated plasma, with a mean free path smaller than an ion gyroradius, v 
is essentially the ion mean free path times an ion thermal speed. It is also the "ion parallel 
viscosity" for the case in which the mean free path is greater than a gyroradius,i@ as it is 
in the case of a tokamak fusion device. For the current generation of tokamaks, v is an ex- 
traordinarily large number, so large as to cast legitimate doubt on the applicability of either 
version of the MHD viscous stress tensor to tokamak dynamics. Considerable discussion 
and some controversy has surrounded the form and magnitude of the appropriate viscous 
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stress tensor to be used in tokamak MHD, and at present, these discussions show no signs of 
converging. We must admit our own doubts about the accuracy of MHD for tokamaks with 
any currently available viscous stress tensor, but for purposes of this discussion, we make 
the assumption that the standard isotropic viscous-stress tensor [used in the derivation of 
Eq. (JJ)] is applicable, but with a viscosity coefficient sufficiently large as to make the viscous 
Lundquist number in Eq. ([!])] small compared to unity: i.e., v is of the order of an 
ion mean free path times an ion thermal speed for the case of a plasma. In any case, the 
assumption made may be satisfied in other magnetofluids and may be considered of possible 
relevance to tokamaks, at least until some more convincing approximation for the viscous 
stress tensor appears. 

This assumption makes it possible to treat the viscous term as of the same order as the 
other two terms on the right hand side of Eq. ([!]), and further makes it possible to neglect 
the inertial terms [the left hand side of Eq. ([[])] for time-independent states. A detailed 
justification in terms of dimensionless numbers is given in the Appendix. The velocity v is 
then to be calculated in terms of j and B from the approximate relation 

Vp-jxB = z/V 2 v, (3) 

where to lowest order, j and B are to be obtained strictly from Ohm's and Ampere's laws 
and the magnetic boundary conditions, without reference to v. The term containing velocity 
in Eq. (|3|) needs to be retained in order to balance the part of the left hand side that has 
a non-vanishing curl. The resistivity and viscosity are to be taken as spatially uniform, so 
that taking the divergence of Eq. (§) would give, using the incompressibility assumption, 

V 2 p = V-(jxB), (4) 

a Poisson equation for the pressure that determines p as a functional of j and B. The 
pressure can also be made to drop out of Eq. (|3]) by taking the curl and writing the result 
as an inhomogeneous equation to be solved for the vorticity uj = VXv: 

Vx (jxB) = -z/V 2 {u v e v ) = -vV 2 u>. (5) 
4 



We shall find that in the geometry considered, the vorticity vector points entirely in the 
(toroidal) ^-direction; see the Appendix. Once the fluid velocity and the vorticity are 
determined by solving Eqs. (|3|) and (0) (v and a> are, as advertised, "small" in the sense of 
being first order in one can return with them to Ohm's and Ampere's laws and iterate 
again, obtaining first-order corrections to the current and magnetic field. Then, corrections 
to v and u) can be obtained by going back again to Eqs. @ and (HD, and iterating. Our 
interest here, though, is in the lowest order solutions for the velocity and vorticity. 

The electric field, highly idealized for tractability, is regarded as being generated by a 
time-proportional axial magnetic field confined to a high permeability cylinder (a cylindri- 
cal iron core, say) whose axis of symmetry is the z-axis, and which extends to infinity in 
the positive and negative z directions. This cylinder lies entirely within the "hole in the 
doughnut" of the toroid and is perpendicular to the mid-plane z = 0. It produces an electric 
field which lies purely in the azimuthal direction [yj-direction in cylindrical polar coordinates 
(r, <f, z), which we use throughout]: E = (E r /r) e v , where E is the strength of the applied 
electric field at a reference radius r = r within the toroid. Ohm's law, in the dimensionless 
units, is (with a resistive Lundquist number 1 / 77) 

E + vxB = r?j, (6) 

so that if we neglect v to lowest order, a purely toroidal current density is generated with the 
form j = (E ro/r]r) e^. The (poloidal) magnetic field associated with this j is determined by 
Ampere's law and by boundary conditions, which we take to be Bn = at the toroidal walls, 
where n is the unit normal. We assume that the toroidal boundaries are highly conducting 
and coated with a thin layer of insulator; we ignore the slits and slots in the conducting 
walls that are required for the applied electric field to penetrate the interior of the toroid: 
a necessary and common if regrettable idealization. We return presently to the explicit 
calculation of the poloidal magnetic field. In addition, we may assume the presence of a 
vacuum toroidal magnetic field that is externally supported: B r = B^e^ = (B r /r) e v . The 
total magnetic field is the sum of the toroidal and poloidal magnetic fields. The toroidal 
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magnetic field seems to play little role in establishing the properties of the equilibrium, 
though it will have a great deal to say about the stability of that equilibrium, a question 
which we do not consider here. The magnetic field lines are, topologically speaking, helical; 
but the toroidal field enters in no other context, and the poloidal magnetic field contains 
all the nontrivial magnetic information. The toroidal magnetic field would also play a much 
more prominent role if a tensor electrical conductivity were allowed. A tensor electrical 
conductivity, which we do not consider, would permit lowest-order poloidal currents as well; 
this limitation would be desirable to remove in future treatments. At the next order, poloidal 
currents would also be implied by the cross product of v with the toroidal magnetic field. 
It is to be stressed that until the velocity field begins to feature in the calculation, the 
mathematics are indistinguishable from a similar low-frequency electrodynamics calculation 
in a finitely conducting ring of metal. It is shown in the Appendix that the condition for 
the negligibility of v in Ohm's law at the lowest order is again the smallness of the viscous 
Lundquist number, not a totally obvious result. 

Before computing the poloidal magnetic field explicitly, we return to Eq. (|5p and examine 
how the vorticity uj may be determined once the left hand side is known. Since the Laplacian 
of a vector in the toroidal direction that only depends on poloidal coordinates is also a vector 
in the toroidal direction that only depends on poloidal coordinates, we essentially have a 
Poisson equation to solve for u> = u^e^. One way to proceed is to expand w in a family of 
vector eigenfunctions of the Laplacian, which are related to waveguide modes: 

V 2 (u v e v ) + A 2 {u) v e v ) = 0. (7) 

All components of all fields are axisymmetric (i.e., the components of all vector fields are 
(/^-independent), and the solution to Eq. (0) is any one of the functions: 



u jk e v = e jk [Ji(jjkr) + D jk Yi(j jk r)] 



i cos kz j 



e 



(8) 



where Ej k , 7^, Dj k and k are undetermined constants, with 7^, k and Xj k related by the 
condition X 2 - k = jj h + k 2 . Here, J\ and Y\ are Bessel and Weber functions, respectively. The 
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vorticity may be written in terms of a velocity stream function ip(r, z) in the following way: 
oJip&ip = — V 2 (ipe^/r). The function ip/r may also be expanded in terms of the eigenfunctions 
obtained from Eqs. (|7|) and The boundary conditions chosen are that ip and u^, and 
thus the ujjk, vanish on the boundary. We now choose a rectangular cross section for the 
toroid; the combined boundary conditions v • n = and u v = make the velocity field 
satisfy stress-free, impenetrable boundary conditions. That is, v • n and the tangential 
viscous stress are zero on all four faces. These are the only boundary-shape and internally 
consistent viscous boundary conditions (perfectly smooth and impenetrable walls) we have 
found for which the solution is obtainable by elementary means. However, we believe the 
qualitative conclusions to be reached will apply to essentially an arbitrary toroidal boundary 
shape. 

We seek a solution to Eq. (|J) by expanding u) v in terms of the eigenfunctions, uijf., having 
determined numerically the allowed values of Ejk, jjk, -Djfc and k. The value of Sjk is chosen 
so that the two-dimensional integral of the square of the eigenfunctions is unity. We assume 
that what we have obtained, then, is a complete orthonormal set. We write 



UJ 



^2fl jk e jk [Jx{jjkr) + D jk Yx{^ jk r)\ 



j,k 



\ cos kz J 



— w i fe ' (^) 

3,k 



with unknown expansion coefficients Qjk- Assuming term-by-term differentiability of the 
expression in Eq. (|D, the stream function is given by 

^ = r ^jk ^jk Ujk- (10) 

j,k 

Using the orthonormal eigenfunctions Ujk, Eq. (5) may be expressed as 

Vx(jxB) = ^A J 2 fc O ifc w jfc e p . (11) 

j,k 

If we choose, we may obtain B from j by writing it in terms of a vector potential B = V X A, 
where A = A v e 9 is to be expanded using the same eigenfunctions as we used to expand the 
vorticity and the stream function. Since the curl of V X A = B yields (if V ■ A = 0) 
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V 2 (Vy = -j v e v = — (12) 

it is clear that A and B can be determined by the previous procedure: solving Eq. ( |12^ ) by 
expanding in the same eigenf unctions (A = on the wall) and differentiating. In the next 
section, we adopt an alternative and slightly simpler procedure to determine B. It is already 
clear, however, what the form is that the vorticity field will have to take, simply by looking 
at the source term in the Poisson equation for vorticity and noting the fact that the radial 
component of B will be positive above the mid-plane and negative below. The vorticity 
distribution will look dipolar in the (r, z)-plane, and the three-dimensional vortices will be 
vortex rings that sit one on top of the other, above and below the mid-plane: a "double 
smoke ring" configuration. 



III. EXPLICIT NUMERICAL SOLUTIONS 

In this section, we determine the poloidal vector fields B p and v for this MHD system. 
These fields are associated with a toroidal current density that is proportional to 1/r, where 
r is the distance from the toroidal axis. For illustrative purposes, we consider a toroid with 
a rectangular cross section (see Fig. 1). We assume the boundaries to lie at the planes 
z = ±L and at the radii r = r_ and r = r + , where r_ < r < r + . In our calculations, 
all length scales are normalized to the major radius ro = (r_ + r+)/2. Since we do not 
assume large aspect ratio, there is no significant difference in the lengths of the major 
and minor radii. Because the walls of the toroid are idealized as perfect conductors, the 
appropriate magnetic boundary condition for this problem is that the normal component of 
the magnetic field vanish there. To find the magnetic field we solve for the magnetic flux 
function x directly. In terms of Xi both the vector potential and the poloidal magnetic field 
can be easily derived. Once B p is known, the source term Vx(jxB) on the left side of 
Eq. (5) can be computed, and the resulting Poisson equation solved for the vorticity. With 
the vorticity, the velocity stream function, ip, and finally the velocity v can also be found. 
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With axisymmetry, the poloidal component of the magnetic field B p may be represented 
in terms of a magnetic flux function x(r, z) according to 

Bp = Vx X V(p. (13) 
Substituting this into Ampere's law, VxB = j, with j = (E r /r]r) e v yields 

A* X = r--^ + ^ = -^. (14) 
dr r dr dz 2 r\ 

Note that the magnetic vector potential A = A^e^ is obtained from the magnetic flux 
function by dividing by r: A v = xjr. We seek x an d B p in the rectangular domain r_ < 
r < r + and — L < z < L subject to the boundary condition B ■ n = 0, where n is the unit 
normal to the wall of the toroid. This boundary condition implies 

= 0, at z = ±L (15a) 

or 
Qx 

— = 0, at r = r + ,r_. (15b) 
oz 



A particular solution of Eq. fli"4}) that vanishes at r = r_ and r = r + is 

= — v ln — - 12 — 12 r + ln — • ( 16 ) 



The solution of the homogeneous equation A*Xh — that is symmetric about the mid-plane 
of the toroid is 

Xh = J2 C K e *> r l J i( Kr ) + («**)] cosh (17) 

re 

where C re , e K , -D K and k are arbitrary constants. The general solution to Eq. (|TJ|) is x — 
X P + Xh- 



Equation ( |15b|) is satisfied by requiring that 



J 1 (Kr_) + D K Y 1 (Kr_) = 0, (18) 
J 1 {Kr + ) + D K Y 1 {Kr + ) = 0, (19) 
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Equations ( |18D and dl9|) can only be solved consistently if the determinant 

V = J 1 (Kr_)Y 1 (Kr + ) - Ji(Kr + )Ki(Kr_) (20) 

vanishes. For an infinite sequence of K-values, with each k corresponding to a particular zero 
of T> for given values of r_ and r + , general Sturm-Liouville theory tells us that the functions 

Ok = e K [J (/cr) + Ac *o(«r)] , (21) 

form a complete ortho normal set on the interval r_ < r < r + . The e K are real constants 
chosen to normalize the 4>q k : 

r+ 

(Pok^ok'T dr = 5 K>K > . (22) 

The z-boundary condition can be satisfied by requiring that 

EfA 1 1 + 2 In — - In T -± \ = *C K e K [ J («r) + A Y (kt)] cosh «L. (23) 

Multiplying both sides of this equation by r 0o/t'( r ) an d integrating from r_ to r + determines 
the coefficients C K . We find 



C. = ^ r + 0o.(r){l + 21n^- 2r + 1 2 n(r+ { r ' ) lrdr/(*coah*L). (24) 
2?] Jr_ { r_ r+ — r: 1 



Table I shows the first ten values of r K, D K) r e K and 2 77 C K / (i?o ?"o) for r_/r = 0.6, r + /r = 
1.4, and L/r = 0.3. With these coefficients and the discrete set of K-values, a magnetic 
flux function x that satisfies both boundary conditions in Eq. (p~5|) can be constructed. A 
plot of the contours of x(r, z) appears in Fig. 2. These contours are the projections of the 
surfaces on which the magnetic field lines lie. Using Eq. (P"3"|), the components of the poloidal 
magnetic field can be easily computed. We find 

B r (r, z) = C K e K n [J t («r) + D K Yi(Kr)] sinh kz, (25) 

K 

B z (r,z) = — — 11 + 2 In ^ In — S 

2?7 [ r - r\ — rt r_ J 

+ ^ C K e K n [Jo{nr) + D K Y (Kr)] cosh kz. (26) 
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The next step in our derivation is to determine the vorticity, uu v . For a current density 
jtp = E r /r]r, the left side of Eq. (11) reduces to — (2E r B r /rir 2 ) e v . Thus, the equation 
to solve is 

E\2 n 2E r B r 
A E \L n £UJ n £ — -, [At) 

where B r is given by Eq. (|25|), and the expansion coefficients fl n e are as of yet undetermined. 
The eigenfunctions u n e that have odd parity in z and vanish on the boundary of the toroid 
are given by 

1 £ttz 

uj ne (r,z) = —= 0i n (r)sin— , £=1,2,3,..., (28) 

y L L 

where \ j\fh is a normalization factor. Possible cos[(2£ + \)kz/2L\ terms may be omitted 
from symmetry considerations. The functions <fii n {r) are defined by 

<Pin(r) = e n [Ji(a n r) + D n Yi(a n r)] , (29) 

where the parameters e n are real constants chosen to normalize 0i n : 

<Pln<Pln> r dr = 5 nin > . (30) 

The "n" in a n designates the n-th zero of the function V, and o? n = \ 2 nl — £ 2 tt 2 /L 2 . Note 
that the a n -values correspond to the K-values, the first ten of which appear in Table I 
for r_/r = 0.6, r + /r = 1.4 and L/r = 0.3. Multiplying Eq. ( [27D by uj u i £ rdrdz and 
integrating over the range r_ < r < r + and — L < z < L determines the expansion 
coefficients Q n e. The result is 

_ 4:7iE r £(-l) e x ^ C m a m sinh a m L 

nl ~ V uL 3 / 2 (a 2 n + £H 2 /L 2 ) ot 2 m + £ 2 k 2 /L 2 

f r + dr 
x / (f>i m (r)<f> ln {r)— . (31) 

A contour plot of ui^ using these coefficients appears in Fig. 3 for r_/r = 0.6, r + /r = 1.4 
and L/r = 0.3. Positive contours are denoted by a solid line and negative contours by a 
dashed line. The vorticity vanishes at the toroidal walls, which is equivalent to stress-free 
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boundary conditions in this problem. The convergence of the vorticity series is rather fast, 
owing to the presence of a 2 + £ 2 n 2 /L 2 terms in the denominator of Q n £. Typically, it was 
only necessary to keep a dozen or so terms to achieve a high degree of accuracy. Contours 
of the velocity stream function ip are shown in Fig. 4 for the same parameters and with the 
same convention. With ip, the velocity can be computed from v = V^xVy. It is easily 
verified that the normal component of the velocity will vanish at the boundaries. In Figs. 3 
and 4, note the appearance of paired- vortex structures that resemble a "double smoke ring" 
configuration. 

Since B r (r, z) in nonzero at z = ±L, a Gibbs phenomenon is to be expected in the series 
for V 2 ^^ near z = ±L. This is a consequence of representing the z-dependence of the 
right side of Eq. ([27]) in terms of sine functions, all of which vanish at z = ±L. Since both 
sides of Eq. ( p7|) are identically zero at r = r_ and r = r + , a Gibbs phenomenon will not 
occur near the boundaries in r. 

IV. DISCUSSION AND CONCLUSIONS 

It is regrettable that nearly all experiments performed during the last several years 
on MHD in toroidal geometry have been carried out in tokamaks intended to confine a 
thermonuclear plasma. Not only are diagnostics for such internal variables as fluid velocity, 
vorticity, and electric current density very limited due to the high temperatures, the detailed 
applicability of MHD itself is in doubt because of previously-mentioned uncertainties in the 
appropriate viscous stress tensor to be used in theoretical models. Here, we have attempted 
to isolate an interesting MHD effect, without taking a position on whether it should or 
should not be an important feature of tokamak operation. Tokamaks are often thought 
to have both overall poloidal and toroidal rotation, which have been attributed to various 
consequences of local charge non-neutrality. The combinations of all three kinds of flows, if 
they were present, might be quite difficult to untangle. 

The "double smoke ring" configuration identified in this paper is a feature associated with 
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electrically-driven toroidal magnetofluids that we believe is quite robust; it does not require 
local charge non- neutrality, and may even appear in liquid metal experiments (for example), 
even though the high-viscosity calculations we have done imply inequalities that may not 
be easily satisfied in liquid metals. Entirely as a consequence of the toroidal geometry, 
a purely toroidal electric current generates a magnetic field for which a part of the j X B 
Lorentz force produces a local toroidally-directed torque on the magnetofluid. (This torque 
disappears in the "straight-cylinder" limit.) This gives rise to opposing pairs of vortex rings 
with vorticity aligned parallel and anti-parallel to the current density. We believe that 
these structures will exist under a variety of boundary conditions (non-conducting walls, 
for example, with no-slip boundary conditions) and will not require low Hartmann numbers 
or viscous Lundquist numbers, though the flows may be more elaborate (involving, say, 
poloidal currents or toroidal velocities as well) when the inequalities we have invoked are 
not satisfied. 

The likely presence of MHD flows in toroidal geometry was probably first reported in an 
unpublished paper by Pfirsch and Schluter! over thirty years ago. Their approach was quite 
different from ours, involving for example an inverse aspect ratio expansion. In addition, 
they ignored the velocity field in the equation of motion (but not in Ohm's law). On the 
basis of their model, Pfirsch and Schluter concluded that there would be a necessary mass 
flux outward from the toroid that required "sources" of mass inside the toroid. Here, we 
have explicitly exhibited a large class of solutions with no normal component of velocity at 
the walls, which contradicts the findings of Pfirsch and Schluter. Nevertheless, credit for 
the observation that flows are to be expected in the steady state must go to them. Ideal 
toroidal vortices have been considered by Marnachev.i 

The flow pattern that we have been computed in this paper, with streamlines that cross 
the toroid near the mid-plane, is not one that is very propitious for plasma confinement 
with high temperatures in the center of the toroid, and lower temperatures near the wall. 
Nevertheless, the ability to "stir" the interior of a toroidal magnetofluid with externally- 
maintained electric fields might have other applications of some interest, such as in the 
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cooling of alloys.1121 More generally, the likely separation of driven MHD states into those 
involving velocity fields and those which are static seems artificial to us. In this particular 
example, we have found velocity fields that do not arise because of any instability or mal- 
function, but are an inherent part of the equilibrium itself, even though no external pressure 
gradients are applied. 

ACKNOWLEDGMENTS 

This work was supported in part by the U.S. Department of Energy under grant DE- 
FG02-85ER53194. 

APPENDIX: REYNOLDS-LIKE NUMBERS AND GEOMETRY 

In this Appendix, we consider in detail the inequalities that justify the neglect to lowest- 
order of v ■ Vv and vxB in Eqs. ([I]) and (|6|), respectively. We also specify the assumed 
toroidal geometry of the fields B, j, v, and uj. 

Proceeding from the dimensional (cgs units) version of Eq. (|T]), the condition for ne- 
glecting the v • Vv term relative to the viscous term is low viscous Reynolds number. This 
Reynolds number is defined as vL/u, where v is a typical fluid velocity, L is a typical length 
scale, and v is a kinematic viscosity. We take the minor toroidal radius as the typical length 
scale for this problem. Using the dimensional version of Eq. (|j|), we may estimate v as C a M, 
where C a is an Alfven speed based on a typical poloidal magnetic-field strength B, and 
where M is the viscous Lundquist number, C a Ljv. Inserting this in the viscous Reynolds 
number requirement, we see that the justification of Eq. (||) follows from the smallness of 
the square of M compared to unity. 

Neglect of the velocity term in Eq. is justified by requiring the velocity v to be small 
compared to rj/L. The parameter rj is the magnetic diffusivity, defined by rj = c 2 /47ra, 
where a is the electrical conductivity. Using the previous estimate for v, we see that this 
inequality is equivalent to SM « 1, where S is the resistive Lundquist number, C a L/rj. 
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Since SM = H 2 , where H is the Hartmann number, an equivalent statement of the second 
inequality is that if 2 is much less than one. Thus, the smallness of the squares of M and H 
compared to unity is enough to justify the approximations made. 

Throughout this paper, we work exclusively in cylindrical polar coordinates (r,(p,z). 
The axis of symmetry is the z-axis, and the mid-plane of the toroid is z = 0. The ip- 
direction is the "toroidal" direction, and the r and z directions are called the "poloidal" 
directions. (See Fig. 1.) For the numerical solutions presented in Sec. Ill, we consider a 
toroid with a rectangular cross section in the (r, z)-plane. The boundaries are idealized 
as perfectly smooth, perfectly conducting rigid walls, with zero normal velocity and zero 
tangential viscous stress. We imagine that the inside surfaces of the boundaries are coated 
with an infinitesimally thin layer of insulator (other idealizations are possible, such as a 
purely non-conducting boundary). In most experimental toroidal devices, gaps are present 
in the conducting boundary to permit the penetration of a toroidal electric field. The 
asymmetric effects introduced by these slits and slots, however, are not included in our 
model. An externally-supported toroidal (vacuum) magnetic field is possible, but will not 
feature in the analysis until a stability calculation or dynamical simulation is performed. 
The nontrivial (poloidal) magnetic field, B p , therefore has only r and z components. The 
current density j has only a (/9-component, the vorticity u) has only a (^-component, and the 
fluid velocity v has only r and z components. It is shown in Sec. Ill that the mechanical 
flow is a (double) "vortex ring" configuration, and the electric current is a (single) "current 
loop." The scalar pressure p will depend only upon r and z. The simple form of the fields in 
this model is a consequence of two assumptions: (i) axisymmetry; and (ii) isotropic scalar 
viscosity and conductivity. The introduction of a tensor conductivity or viscosity moves the 
problem out of the realm of present tractability. 

Note that the part of the j X B term in Eq. (|3]) that contains a curl is of first order in the 
inverse aspect ratio (i.e., the ratio of minor and major radii of the toroidal system). The 
results of this paper, however, are not predicated on the assumption that the inverse aspect 
ratio is small compared to one. 
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FIG 1 . Geometry of computational model. The toroid has 
a rectangular cross section with impenetrable, perfectly 
conducting, perfectly smooth walls. The magnetofluid 
occupies the region between the radii r_ and r + , and the 

planes z = ± L. 




FIG. 2. Contours of the magnetic flux function %(r,z) for 
r_ / r Q = 0.6, r + / r Q = 1 .4 and L / r Q = 0.3. 
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FIG. 3. Contours of the toroidal vorticity co^nz) for 

r_ / r Q =0.6, r + / r Q =1 .4 and L / r Q = 0.3. Positive 

contours are denoted by a solid line and negative 
contours by a dashed line. 
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FIG 4. Contours of the velocity stream function \|/(r,z) for 
r_ / Tq = 0.6, r + / Tq = 1 .4 and L / rg = 0.3. Positive contours 

are denoted by a solid line and negative contours by a 
dashed line. 



